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Abstract 

Most extensions of the Standard Model lepton sector predict large lepton flavor 
violating rates. Given the promising experimental perspectives for lepton flavor 
violation in the next few years, this generic expectation might offer a powerful 
indirect probe to look for new physics. In this review we will cover several 
aspects of lepton flavor violation in supersymmetric models beyond the Minimal 
Supersymmetric Standard Model. In particular, we will concentrate on three 
different scenarios: high-scale and low-scale seesaw models as well as models 
with R-parity violation. We will see that in some cases the LFV phenomenology 
can have characteristic features for specihe scenarios, implying that dedicated 
studies must be performed in order to correctly understand the phenomenology 
in non-minimal supersymmetric models. 


1 



Contents 


1 Introduction 3 

2 Current experimental situation and future projects 5 

3 High-scale seesaw models 6 

3.1 Standard high-scale seesaw scenarios. 8 

3.2 Extended high-scale seesaw scenarios. 11 

4 Low-scale seesaw models 15 

4.1 The supersymmetric inverse seesaw. 16 

4.2 LEV in the supersymmetric inverse seesaw. 17 

5 R-parity violating models 19 

5.1 Explicit R-parity violation. 21 

5.2 Spontaneous R-parity violation . 25 

6 Summary and conclusions 27 


2 








1 Introduction 


The Standard Model (SM) particle content has been recently completed with the discovery 
of the long-awaited Higgs boson at the CERN Large Hadron Collider (LHC) [1,2]. This 
constitntes a well deserved reward after decades of intense search, with great efforts from 
the theory and experimental commnnities. Fnrthermore, it also confirms that the SM must 
be, at least to a good approximation, a precise description of nature up to the energies 
explored. In fact, and apart from some phenomenological facts that indeed require some 
unknown new physics (NP), like the existence of dark matter and neutrino masses, the 
SM explains to a high level of accuracy all the observations made in a wide variety of 
experiments. 

For the last decades, the progress in theoretical particle physics has been driven by 
naturalness considerations in the form of the famous hierarchy problem. This has led to 
many extensions of the SM, all of them attempting to explain why the weak scale has not 
been pushed to much higher energy scales by some hypothetical NP degrees of freedom. 
Among the many proposals to address this issue, supersymmetry (SUSY) is certainly the 
most popular one. However, and similarly to other analogous solutions to the hierarchy 
problem, the predicted new particles at the weak scale have not been observed at the LHC. 

This has of course raised some doubts about the existence of supersymmetry close to the 
weak scale. Since this proximity is to be expected in case supersymmetry has something to 
do with the hierarchy problem, the whole idea of weak scale supersymmetry is under some 
pressure at the moment. However, it is worth keeping in mind that most experimental 
searches for SUSY focus on the Minimal Supersymmetric Standard Model (MSSM). This 
model, which constitutes the minimal extension of the SM that incorporates SUSY, has 
some underlying assumptions that lead to very specific signatures. For example, in the 
MSSM one assumes the conservation of a discrete symmetry, known as R-parity [3,4], 
which forbids all renormalizable lepton and baryon number violating operators and leads 
to the existence of a stable particle, the lightest supersymmetric particle (LSP), which 
in turn leads to large amounts of missing energy in supersymmetric events at the LHC. 
There are, however, many known (and well motivated) supersymmetric scenarios with 
R-parity violation and, in fact, several authors have shown that simply by allowing for 
non-zero B-violating terms in the superpotential, the current LHC bounds can be clearly 
relaxed, allowing for the existence of light squarks and gluinos hidden in the huge QCD 
background [5,6]. Similarly, one can extend the MSSM in many other directions, often 
changing the phenomenology at colliders dramatically. This suggests that it might be too 
soon to give up on SUSY, a framework with many possibilities yet to be fully explored. 

As explained above, there are some well-grounded phenomenological issues that cannot 
be explained within the SM. One of these open problems is the existence of non-zero 
neutrino masses and mixings, nowadays firmly by neutrino oscillation experiments [7-9]. 
In fact, this issue is not addressed in the MSSM either, since neutrinos remain massless 
in the same way as in the SM. This calls for an extension of the MSSM that extends the 
lepton sector and accommodates the observations in neutrino oscillation experiments. This 
can be done in two different ways: (1) high-energy extensions, in which the new degrees of 
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freedom responsible for the generation of neutrino masses live at very high energy scales, 
and (2) low-energy extensions, with new particles and/or interactions at the SUSY scale. 

One of the most generic predictions in neutrino mass models is lepton flavor violation 
(LFV). In fact, neutrino oscillations are the proof that lepton flavor is not a conserved sym¬ 
metry of nature, since neutrinos produced with a given flavor change it as they propagate. 
Therefore, all neutrino mass models built to give an explanation to oscillation experiments 
violate lepton flavor. However, we have never observed LFV processes involving charged 
leptons although, in principle, there is no symmetry (besides lepton flavor, which we know 
to be broken) that forbids processes like fi~ —)■ e“ 7 , t~ —)■ or —)■ e~fi~^. This 

fact can be well understood in some minimal frameworks, such as the minimal extension 
of the SM with Dirac neutrinos. In this case, LFV in the charged lepton sector is strongly 
suppressed, since neutrino masses are the only source of LFV, leading to unobservable 
LFV rates, like BR(/i“ —)■ e~'y) ~ 10“^® [10]. However, as soon as one extends the SM, 
this conclusion can be clearly altered [11,12]. In fact, new sources of LFV can be found in 
most extensions of the leptonic sector, either caused by new interactions, by new particles 
or even by complete new sectors that couple to the SM leptons. 

After this discussion on LFV and neutrino masses a clarihcation is in order. Although 
neutrino oscillations imply LFV, LFV does not necessarily imply neutrino oscillations. 
There are models that predict charged lepton LFV without generating a mass for the 
neutrinos. The simplest example of this class of models is the general Two-Higgs-Doublet of 
type-HI, where neutrinos remain massless but lepton flavor is violated due to the existence 
of off-diagonal h — ii — ij vertices. Another relevant example is the MSSM itself, where 
neutrinos are also massless, but the slepton soft masses can induce LFV processes if they 
contain off-diagonal entries. One can actually estimate the branching ratio for the radiative 
LFV decay ii —>■ ij'y as [13] 


BR(£, ^ £,7) ^ U BR(£, ^ i^UiPj ), (1) 

^^^SUSY 

where Gp is the Fermi constant, a the hne structure constant, are the dominant 

off-diagonal elements of the soft SUSY breaking slepton mass matrices and Msusy is the 
typical mass of the SUSY particles, expected to be in the TeV ballpark. This estimate 
clearly shows that rather small off-diagonal elements are required to satisfy the experimen¬ 
tal bounds [14]. 

In general, one concludes that large LFV rates are expected in most models beyond the 
SM. This motivates the study of LFV as an indirect probe of new physics and, in particular, 
of supersymmetric models beyond the MSSM. This is the subject of this review In 
particular, we will concentrate on three different scenarios: high-scale and low-scale seesaw 
models as well as models with R-parity violation. As we will see, the LFV phenomenology 
turns out to be very different depending on the exact scenario, implying that lepton flavor 

^The field of lepton flavor violation beyond the MSSM has been intensely explored for many years 
and contains a vast literature. In this review I present my personal view of the subject and thus I must 
apologize for those papers which are not cited. 
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violation may be richer than in the MSSM. In some cases the common lore (established 
in the MSSM) turns out to be wrong, and specihc studies must be performed in order to 
correctly understand the corresponding LFV phenomenology. 

This review is organized as follows: in Sec. 2 we give an overview of the current exper¬ 
imental situation and briefly discuss some projects that will take place in the near future. 
Then we review the LFV phenomenology of three different types of models beyond the 
MSSM: high-scale seesaw models (in Sec. 3), low-scale seesaw models (in Sec. 4) and 
models with R-parity violation (in Sec. 5). Finally, we conclude in Sec. 6. 


2 Current experimental situation and future projects 

The search for LFV is soon going to live a golden age given the upcoming experiments 
devoted to high-intensity physics In addition to the LFV searches already taking place 
in several experiments, new projects will join the effort in the next few years. 

In what concerns the radiative decay ^i —)■ ij'y, the experiment leading to the most 
stringent constraints is MEG. This experiment, located at the Paul Scherrer Institute in 
Switzerland, searches for the radiative process /i —)■ ey. Recently, the MEG collaboration 
announced a new limit on the rate for this process based on the analysis of a dataset 
with 3.6 X 10^^ stopped muons. The non-observation of the LFV process led to the limit 
BR(/i —)■ ey) < 5.7 ■ 10“^^ [18], four times more stringent than the previous limit obtained 
by the same collaboration. Moreover, the MEG collaboration has announced plans for 
future upgrades. These will allow to reach a sensitivity of about 6 ■ after 3 years of 

acquisition time [19]. This is of great importance, as this observable is the one with the 
largest rates in many models. 

The most promising improvements in the near future are expected in p —)• 3 e and fi — e 
conversion in nuclei. Regarding the former, the decay /i —)■ 3 e was searched for long ago 
by the SINDRUM experiment [20], setting the strong limit BR(/i —)■ 3e) < 1.0-10“^^. The 
future Mu3e experiment announces a sensitivity of ~ 10“^® [21], which would imply an 
impressive improvement by 4 orders of magnitude. As for the latter, several experiments 
will compete in the next few years, with sensitivities for the conversion rate ranging from 
10“^^ to an impressive 10“^®. These include Mu2e [22-24], DeeMe [25], GOMET [26,27] 
and the future PRISM/PRIME [28]. In all cases, these experiments will dehnitely improve 
on previous experimental limits. 

The limits for r observables are less stringent, although signihcant improvements are 
expected from B factories like Belle II [29,30]. Finally, although the most common way 
to search for LFV is in low-energy experiments, colliders can also play a very relevant role 
looking for LFV processes at high energies. The LHGb collaboration reported recently 
the hrst bounds on r —?• 3p ever obtained in a hadron collider [31]. Furthermore, the 
GMS collaboration recently found an intriguing 2.4(T excess in the h —?■ rp channel which 
translates into BR(h —)■ r/i) = (0.841 q; 37)% [32]. For reference, in Tab. 1 we collect 

^See [15-17] for recent reviews. 
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LFV Process 

Present Bound 

Future Sensitivity 

p —)■ ey 

5.7 X 10-^3 [ig] 

6 X 10-1^ [19] 

r —)■ ey 

3.3 X 10-® [33] 

~ 3 X 10-9 [29] 

r —)■ /iy 

4.4 X 10“® [33] 

~ 3 X 10-9 [29] 

fi — >■ eee 

1.0 X 10-^2 [20] 

~ 10-1® [21] 

T —>■ /i/i/i 

2.1 X 10“® [34] 

~ 10-9 [29] 

T~ —>■ 

2.7 X 10-® [34] 

~ 10-9 [29] 

T~ —)■ p“e+e“ 

1.8 X 10-® [34] 

~ 10-9 [29] 

r —)■ eee 

2.7 X 10-® [34] 

~ 10-9 [29] 

Ti — )■ e“, Ti 

4.3 X 10-12 [35] 

~ 10-13 [36] 

/i“, Au —)■ e~ , Au 
p , A1 —^ e , A1 
/i-,SiC ^ e-,SiC 

7 X 10-13 [37] 

10-1® - 10-13 
10-1^ [38] 


Table 1: Current experimental bounds and future sensitivities for the most important LFV 
observables. 

present bounds and expected near-future sensitivities for the most popular low-energy 
LFV observables. 

The theoretical understanding of all these processes will be crucial in case a discov¬ 
ery is made. With such a large variety of processes, the determination of hierarchies or 
correlations in specihc models will allow us to extract fundamental information on the un¬ 
derlying physics behind LFV. This goal requires detailed analytical and numerical studies 
of the different contributions to the LFV processes, in order to get a global picture of the 
LFV anatomy of the relevant models and be able to discriminate among them by means 
of combinations of observables with dehnite predictions [39]. 


3 High-scale seesaw models 


Neutrino mixing is, by itself, a flavor violating effect. Therefore, all neutrino mass models 
that aim at explaining the observed pattern of neutrino masses and mixings incorporate 
lepton flavor violation. However, specihc predictions can be very different in different 
models. 

Among the huge number of scenarios proposed for neutrino mass generation, the see¬ 
saw mechanism is arguably the most popular one. In its conventional form, the seesaw 
mechanism explains the smallness of neutrino mass by means of a very large energy scale, 
the seesaw scale Mss, which suppresses neutrino masses as 


Mss 


( 2 ) 


Here (iL°) = v/\/2 = 174 GeV is the standard Higgs boson vacuum expectation value 
(VEV) that determines the weak scale. In order to obtain neutrino masses of about ~ 0.1 
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eV, one requires Mss ~ 10^^ GeV. For this reason, this setup is usually called high-scale 
seesaw. The proximity of the high-energy scale Mss to the grand unihcation (GUT) scale 
(as predicted in the MSSM) tugut = 2 ■ 10^® GeV suggests an intriguing connection with 
unihcation physics, making the seesaw a very well-motivated scenario. 

Regarding specihc realizations of the seesaw mechanism, it is well-known that with 
renormalizable interactions only, three tree-level realizations exist [40]. These are usually 
called type-I [41-46], II [45-51] and III [52]. They differ from each other by the nature of 
the seesaw messengers: in the type-I seesaw these are singlet right-handed neutrinos, in 
the type-II seesaw scalar SU{2)l triplets with hypercharge two, and in the type-III seesaw 
fermionic SU{2)l triplets with vanishing hypercharge. In all cases they lead to a neutrino 
mass of the form of Eq. (2), where Mss is proportional to the mass of the heavy mediators, 
and the induced neutrino masses are of Majorana type, thus breaking lepton number in 
two units. 

Given the large Majorana masses of the seesaw mediators, one may wonder about how 
to probe the high-scale seesaw. In supersymmetric scenarios this is possible thanks to 
the sleptons. Even if their soft terms are flavor conserving at some high-energy scale, the 
renormalization group running down to the SUSY scale will induce non-zero off-diagonal 
terms due to their interactions with the seesaw mediators [53]. These can be probed since 
the misalignment of the slepton mass matrices with respect to that of the SM charged 
leptons induces LEV processes such as U —)■ ijj, U —>■ 3£j and /i — e conversion in nuclei. 
This connection between the phenomenology at low-energies and the high-scale mediators 
is only possible in supersymmetric models and constitutes an excelent opportunity to test 
the standard seesaw scenario 

In case of high-scale seesaw models, the low-energy theory is simply the MSSM. This 
allows one to establish dehnite patterns and hierarchies among the LEV observables. For 
instance, the branching ratios for the U —>■ ij'y and U —t 3ij LEV decays follow the 
approximate relation [55-57], 



( 3 ) 


Therefore, in supersymmetric high-scale seesaw models, the most constraining LEV process 
is ii —)■ ij'y. The relation in Eq. (3) is caused by the so-called dipole dominance in high- 
scale seesaw models. Among the different contributions to the 3-body decay U —)■ 3£j, the 
dipole photon penguins are the dominant ones, leading to the proportionality between the 
G —t £jy and —)■ 3£j branching ratios 

^In the non-SUSY version of the seesaw mechanism this link between high and low energy scales is lost. 
In this case probing the origin of neutrino masses becomes a quite challenging task, and only very indirect 
probes such as neutrinoless double beta decay are possible [54]. 

An exception to this general rule is found for low pseudoscalar masses and large tan [58] . In this 
case, Higgs penguins turn out to be dominant in processes involving the second and third generations, like 
T —> 3 /r. However, this region of parameter space is nowadays under some tension due to strong flavor 
constraints derived from the observation of quark flavor violating processes like Bg —>■ gd' p~ [59]. 
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3.1 Standard high-scale seesaw scenarios 

Implementing a high-scale seesaw mechanism in supersymmetric scenarios involves an ad¬ 
ditional complication. This is related to one of the most appealing features of the MSSM: 
gauge coupling unihcation. In case of the type-I seesaw, the introduction of the seesaw 
mediator does not spoil this attractive feature, since the right-handed neutrino superhelds 
are gauge singlets and do not affect the running of the gauge couplings. In contrast, in the 
type-II and type-III seesaws, new contributions to the running of the SU{2)i and U{1)y 
gauge couplings are induced by the seesaw mediators. However, a well-known solution to 
this problem exists. Unihcation can be easily restored by embedding the seesaw mediators 
in full SU{5) multiplets, like 15-plets in the case of type-II [60] or 24-plets [61] in the 
case of type-III. The contributions from the other members of the multiplet guarantee that 
the three gauge couplings will eventually meet at a high energy scale, mcuT, although the 
common value of the coupling changes, Qgvt, might be different from that of the MSSM. 
In addition, note that the 24-pIet of SU{5) contains, besides the SU{2)l triplet, a singlet 
state which also contributes to neutrino masses. Hence, in this case one actually has a 
mixture between type-I and type-III seesaws. 

The new superheld content, explicitly denoting gauge charges under SU (3)c x SU (2 )l x 
U{1)y, and superpotential for each seesaw variant are [62]: 

• Type-I : Three generations of right-handed neutrino superhelds, singlets of SU{5), 

are introduced, (1,1,0). 

Wj = IUmssm + (4) 

• Type-II: In this case one needs to introduce a vector-like pair of 15 and 15 of 
S'f/(5), decomposed as S' ~ (6,1, —2/3), T ~ (1, 3,1) and Z ~ (3, 2,1/6) (as well as 

the corresponding bar superhelds). T and T are the SU{2)l triplets responsible for 
neutrino mass generation. 


Wii = Wmssm +^{YTLfL + Ysd‘^Sd‘^) + Yzd‘^ZL 

v2 

+ -^{X,HdfHd + X2HuTHu)+MTfT + MzZ^ + MsS^ (5) 

v2 

• Type-III: Three generations of 24 of SU{5) are added. They can be decomposed as 

~ (1,1, 0), G ~ (8,1, 0), E ~ (1,3, 0), X ~ (3, 2, -5/6) and X ~ (3, 2, 5/6). As 
explained above, neutrino masses are generated as a combination of a type-I seesaw 
(mediated by X'^) and a type-III seesaw (mediate by the SU{2)l triplet E). 



Figure 1: Branching ratios for ^i —>■ £^7 and li —>■ ?>^j as a function of the seesaw scale 
in the SUSY type-I seesaw. This hgure was obtained in the standard SPSla’ point, as¬ 
suming degenerate right-handed neutrinos and hxing the neutrino Yukawas to reproduce 
tribimaximal mixing. Furthermore, a massless lightest neutrino was also assumed. Figure 
taken from [70]. 


Wui = IUmssm + Hu I IeS - \l-YuN^ ) L + YxE^Xd^ 


+ Ylf,N^N‘’ + YicGG + iAfsEE + MxXX 


( 6 ) 


The following notation is used in Eqs. (4), (5) and ( 6 ): IFmssm is the MSSM superpo¬ 
tential, Hd, Hu and L are the down-Higgs, up-Higgs and lepton SU{2)l doublet superfields, 
respectively, and (f is the right-handed down-type quark superfield. 

The LFV phenomenology of SUSY seesaw models has been studied by many authors. 
For the type-I seesaw, low-energy LFV decays such as U —)■ f'j-y and U —>■ Sf'j have been 
calculated in [56,57,63-73]. Similarly, /i —e conversion in nuclei has been studied in [74,75]. 
The other two seesaw variants have received much less attention. The LFV phenomenology 
of the SUSY type-II seesaw has been considered in [60,76-81], whereas the SUSY type-III 
seesaw has been studied in [82-84]. More recently, the interplay between the Higgs mass 
constraint and LFV was studied in [62] for the three seesaw variants. In the following we 
comment on some selected results. 

Let us hrst comment on some results for the SUSY type-I seesaw. Fig. 1 shows the 
ranching ratios for the li —)■ £j 7 and —)■ 3£j decays as a function of the seesaw scale 
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(the mass of the right-handed neutrino mass). This hgure was obtained in Ref. [70], 
using the standard SPSla’ point [85], assuming degenerate right-handed neutrinos and 
a massless lightest neutrino and fixing the neutrino Yukawas to reproduce tribimaximal 
mixing. Although this parameter choice is nowadays excluded for several reasons (the 
SUSY spectrum is too light to pass the constraints from LHC searches and tribimaximal 
mixing is now excluded after 6^13 has been measured), it serves to illustrate the dipole 
dominance discussed above. Indeed, one sees a perfect correlation between the branching 
ratios of ii —)■ ij'y and ii —)■ 3ij, with BR(£j —)■ 3ij) BR(C —)■ 3ij). As already 
discussed, this is due to the fact that the photonic dipole operator dominates 

both processes. 

We now turn ot the SUSY type-II seesaw. In the type-II seesaw, the neutrino mass 
matrix is proportional to the Yt Yukawa matrix. 




2 Mx 


Yt. 


(7) 


This is derived from the superpotential term YtLTL in Eq. (5). This direct relation has 
important consequences for the phenomenology, since it forces the flavor structure of Yt to 
be the same as that of the latter being measured in neutrino oscillation experiments 
In other words, if all the neutrino masses, angles and phases were known, Yt would be 
completely fixed (up to an overall constant). Since Yt determines the LEV phenomenology, 
this implies correlations between the neutrino oscillation parameters and LEV observables. 

A clear illustration of the previous point is shown in Fig. 2, borrowed from Ref. [77]. 
By computing the ratios BR(£j —)■ ij'y) /BR(f'm —t C 7 ) one gets rid of the unknown overall 
factor in the Yt Yukawas, thus obtaining direct predictions in terms of neutrino parameters. 
In this case, the hgure shows the dependence of these ratios on the mixing angle ^13 and 
the Dirac CP violating phase 6. We see that this scenario is extremely predictive. For 
example. Ending experimentally BR(r —)■ ey) > BR(/i —)■ ey) would immediately rule out 
the model, at least in its minimal form 

Additional ways to test high-scale SUSY seesaws include slepton mass splittings [87] 
(directly related to LEV) and the study of the SUSY spectrum, usually deformed with 
respect to the standard spectra in constrained (CMSSM) scenarios. In particular, one can 
construct certain invariants that contain information about the high-energy scale, see for 
example [61,88,89]. See also [90] for related ideas. 

®In contrast, in the type-I and type-III seesaws the analogous relation is quadratic in the Yukawa 
coupling. This introduces extra freedom in the determination of the seesaw parameters (usually encoded 
in the so-called R matrix [86]) and makes it impossible to predict the Yukawa flavor structure only from 
neutrino oscillation data. 

®One way to spoil these strict predictions is to introduce a second SU{2)l triplet T'. In this case m^, 
would receive contributions from T and T', = mf + mf , and the proportionality in Eq. (7) would be 

lost. 
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BR.-r-BR-r^c"; /BR,/i->e7 



■513 


BRr—>^7/BR,p_>f.- BR.-r^e7/BRyt—>e7 



•513 


Figure 2: Contours of the ratios BR(r —)■ /i7)/BR(/i —)■ 67) (black, dashed lines) and 
BR(r —)■ e7)/BR(/i —)■ 67) (red, solid lines) in the (sin 6 * 13 , 5) plane, for normal hierarchy 
(on the left) and inverted hierarchy (on the right) for the neutrino mass spectrum. Figure 
taken from [77]. 


3.2 Extended high-scale seesaw scenarios 


We now turn our attention to extended high-scale SUSY seesaw scenarios beyond the 
classical type-I, type-II and type-III seesaws. However, before we concentrate on the 
extended models, let us make a general observation. As already discussed, flavor violating 
entries in the slepton soft terms m| and m? (the left and right slepton squared soft masses, 
respectively) are induced due to their interactions with the seesaw mediators. Even if they 
are flavor diagonal at the unihcation scale, off-diagonal terms are generated at low energies 
by renormalization group running, thus inducing all kinds of LFV processes. In the case 
of the radiative ii —)■ ij'y, the effective dipole operator that contributes to the decay can 
be written as 

-Cdipoie = e ^ + h.c., (8) 

where Pl,r = |( 1 Y 75 ) are the usual chirality projectors and e the electric charge. The Wil¬ 
son coefficients and Aji are generated by loops with left and right sleptons, respectively. 
One hnds 


Al¬ 


im 




MIusy 






^^^SUSY 


(9) 


where it has been assumed that A-terms mixing left-right transitions are negligible. BR(U —>■ 
£j 7 ) can be computed in terms of and Ar as 


BR(U —)■ £j 7 ) ~ „ 2 — fl^/1 I BR(U —>■ 1 


Gl 


( 10 ) 


The straightforward combination of Eqs. (9) and (10) leads to Eq. (1). 
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Superheld 

generations 

SU{3% 

SU{2)l 

SU{2)r 

U{1)b-l 

Q 

3 

3 

2 

1 

1 

3 


3 

3 

1 

2 

1 

3 

L 

3 

1 

2 

1 

-1 


3 

1 

1 

2 

1 

$ 

2 

1 

2 

2 

0 

A 

1 

1 

3 

1 

2 

A 

1 

1 

3 

1 

-2 

A'’ 

1 

1 

1 

3 

-2 

A 

1 

1 

1 

3 

2 

h 

1 

1 

3 

1 

0 


1 

1 

1 

3 

0 


Table 2: LR model. Matter content between the GUT scale and the SU{2)pi breaking 
scale. The electric charge operator is dehned as Q = hi + hn + 

In the minimal SUSY seesaw models discussed above, the seesaw mediators only couple 
to the left sleptons. For instance, in the type-I case this interaction is given by the super¬ 
potential coupling YuLH^N , whereas in the type-II case it is given by the Y^LTL term. 
For this reason, negligible off-diagonal entries in m| are induced, implying that minimal 
SUSY seesaw models predict Ar ~ 0. As we will see below, this has an impact on some 
low-energy observables that allow, in principle, to test the minimality of the high-scale 
seesaw mechanism. 

Supersymmetric models with non-minimal seesaw mechanisms 

As an example supersymmetric model with a non-minimal seesaw mechanisms, we consider 
the left-right symmetric model of [91,92] (in the following simply called ‘the LR model’). 
The LFV and dark matter phenomenology of this model has been studied in detail in 
[93,94]. 

The model is dehned below the GUT scale where the gauge group is SU{3)c x 
SU (2)£ X SU{2)ji X U (1) b-l- In addition, we assume that parity is conserved. The matter 
content of the model is given in Tab. 2. Here Q, Q^, L and L'’ are the quark and lepton 
superhelds of the MSSM with the addition of (three) right-handed neutrino superhelds to 
complete the L'’ SU{2)b doublets. 

Two <F superhelds, bidoublets under SU{2)l x S'U(2)r, are introduced. Among their 
components, they contain the standard and Hu MSSM Higgs doublets. Finally, the 
rest of the superhelds in Tab. 2 are introduced to break the LR symmetry. 

^The model implicitly assumes the existence of a GUT model at higher energies. At mcuT, the gauge 
couplings and soft terms unify. 
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With the representations in Tab. 2, the most general snperpotential compatible with 
the gange symmetry and parity is 


Wlr = YQQ^Q’^ + YLL^L^-f^M + fLAL + fL^A^L^ 

+ aAhA + 

+ MaAA + M^A'^AV (11) 


Family and gange indices have been omitted in Eq. (11), more detailed expressions can be 
fonnd in [91]. Note that this snperpotential is invariant under the parity transformations 
Q ^ {Qy, L o {Zy, $ ^ A ^ (Ay, A ^ (A)*, Q o (ll")*. This discrete 
symmetry reduces the number of free parameters of the model. 

The breaking of the left-right gauge group to the MSSM gauge group takes place in 
two steps: SU(2)ji x U(1)b-l — t U(1)r x U(1)b-l — t U(1)y- In the hrst step, the neutral 
component of the triplet hi takes a VEV, 

{SJ'») = ^, (12) 

which breaks SU{2)b- However, since = 0 there is a U{1)b, symmetry left over. 

Next, the group U(1)r x U(1)b-l is broken by 


(A"°) 


vbl 

V2 ’ 


(A"°) 


Vbl 

V2 ■ 


(13) 


The remaining symmetry is now U(1)y with hypercharge dehned as E = hn + ^^2^- 

Regarding neutrino masses, assuming that the left triplets (A and A) have vanishing 
VEVs, one induces neutrinos masses from a type-I seesaw only thanks to the presence of 
the right-handed neutrinos [91]. 

Before discussing how to test this scenario with lepton flavor violation, let us mention 
some other non-minimal SUSY seesaw models. The phenomenological study in [95] is 
based on a model very similar to the discussed here, without H superhelds. See also 
[96] for a comprehensive study of supersymmetric models with extended gauge groups at 
intermediate steps. Finally, the seesaw mechanism can also be embedded in SUSY GUTs, 
usually leading to very predictive scenarios [97-102]. 


Probing non-minimal seesaw mechanisms 

As already discussed, a pure seesaw model predicts Ar — 0 simply because the right 
sleptons do not couple to the seesaw mediators. However, in models with non-minimal 
seesaw mechanisms, new interactions between the right sleptons and the members of the 
extended particle content at high energies might exist. When this is the case, non-zero Ar 
coefficients can be induced. 

Let us consider an example. In the LR model, the left-right symmetry implies that, 
above the parity breaking scale, the flavor violating entries generated in are exactly as 
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Figure 3: Positron polarization asymmetry ^(/i+ —>■ as a function of the ratio -Rblr = 

log(uij/mGUT)/log(usL/mGUT)- The seesaw scale Mss has been hxed to 10^^ GeV, whereas 
vbl and vr take values in the ranges vbl ^ [10^"^, 10^^] GeV and vr G [10^^, 10^®] GeV. 
Lighter colours indicate larger vrl- The GMSSM-like parameters have been taken as in 
the SPSS benchmark point [85]. Figure taken from [93]. 


large as the ones in m|. As a consequence of this, Ar 7 ^ 0 is obtained at low energies. In 
fact, one can even get a handle on the symmetry breaking pattern at high energies. Below 
the SU{2)r breaking scale, parity is broken and left and right slepton soft masses evolve 
differently. The left ones keep running from the SU{2)r breaking scale to the U{1)b-l 
scale due to the left slepton couplings with the right-handed neutrinos. One thus expects 
larger flavor violating effects in the left slepton sector, and the difference between left and 
right must correlate with the ratio vrl/vr, which measures the hierarchy between the two 
breaking scales. 

The question is how to measure this difference. For this purpose one can use the 
positron polarization asymmetry, defined as 




\Al\‘^-\Ar\^ 
\As\^ + \Ar\^' 


(14) 


If MEG observes /i"*" —>■ e’''y events, the angular distribution of the outgoing positrons 
can be used to discriminate between left- and right-handed polarized states and measure 
A [103,104]. And this can in turn be used to get information on Ar and Ar. 

In a pure SUSY seesaw model one expects M ~ -|-1 to a very good accuracy. However, 
in models with non-minimal seesaw mechanisms A can signihcantly depart from -|- 1 . For 
example, the LR model typically leads to significant departures from this expectation, 
giving an interesting signature of the high-energy restoration of parity. This is shown 
in Fig. 3, extracted from [93]. First of all, it is clear than the polarization asymmetry 
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—)■ e+ 7 ) is well correlated with the quantity log(ui?/mGUT)/log(uBL/''^GUT)- One 
hnds that as vbl and vr become very different, A approaches +1. In contrast, when 
the two breaking scales are close, vbl/vr ~ 1 , this effect disappears and the positron 
polarization asymmetry approaches ^ = 0 . Note that a negative value for A is not possible 
in this model, since the LFV terms in the right slepton sector never run more than the 
corresponding terms in the left one. 

There are alternative ways to test non-minimal high-scale SUSY seesaw scenarios. 
These include the study of the SUSY spectrum and, in particular, of the invariants pointed 
out for minimal seesaw models. In this case, they contain information about the high- 
energy intermediate scales [96,105] 

4 Low-scale seesaw models 

The high-scale seesaw has an important drawback: the heaviness of the seesaw mediators 
precludes any chance of direct tests. Only indirect tests, based on low-energy processes 
which may have an imprint of the high seesaw scale Mrs, are possible, as explained in 
Sec. 3. In contrast to high-scale models, low-scale seesaw models [106] offer a richer 
phenomenological perspective since the seesaw mediators are allowed to be light. In this 
type of neutrino mass models, instead of Eq. (2), neutrino masses are given by 

~ A. ^ . (15) 

where Mrs is again given by the mass scale of the seesaw mediators and (not to be 
confused with the /i parameter of the MSSM) is a small dimensionful parameter, <C 
V, Mrr. In this case, the smallness of neutrino masses is not obtained with a large Mrr 
scale, but with a tiny parameter. Indeed, if Mrr ~ TeV, Eq. (15) implies a 
parameter of the order of the eV in order to get ruj, in the ~ 0.1 eV ballpark. Therefore, 
one can simultaneously obtain the correct size of neutrino masses while having seesaw 
mediators at the TeV scale. This leads to a plethora of new effects, not present in high- 
scale seesaw models, induced by the light seesaw mediators. In particular, novel (and 
sizable) contributions to LFV processes are possible, sometimes breaking the relation in 
Eq. (3). 

The /ij, parameter is intimately related to the breaking of lepton number. In fact, in 
the /ij, —)■ 0 limit, lepton number is restored and the Majorana neutrino masses in Eq. (15) 
vanish. This makes the smallness of the parameter natural, in the sense of’t Hooft [107], 
since the symmetry of the Lagrangian gets increased when the parameter is set to zero. 
For this reason, low-scale seesaw models are also said to have almost conserved or slightly 
broken lepton number. 

The collider phenomenology of low-scale seesaw models is much richer than that of high- 
scale ones. The seesaw mediators can in principle be produced and, through their decays, 
one may be able to test the mechanism behind neutrino masses. At the LHC, one typically 
expects multilepton final states, often including missing energy carried away by undetected 
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neutrinos. In addition, the LFV signatures can be as frequent as the flavor conserving 
ones. For an imcomplete list of references on the phenomenology of low-scale seesaw 
models see [108-131]. In the case of a type-I seesaw, the seesaw mediator is a fermionic 
gauge singlet. This usually suppresses its production in hadronic colliders. However, 
sizable right-handed neutrino production cross-sections are possible in some type-I seesaw 
realizations due to the mixing with the left-handed neutrinos, which serves as a portal 
to the gauge sector. Furthermore, when the type-I seesaw is embedded in a left-right 
symmetric scenario [132-134] new production mechanisms are possible thanks to the new 
charged currents mediated by the gauge bosons. This allows for further collider tests 
or the model, including searches for lepton number violation, see for example [135-141]. 

We now present the most popular representative of the low-scale seesaw models: the 
inverse seesaw. For other low-scale seesaw models and their LFV phenomenology see 
[142-147]. 

4.1 The supersymmetric inverse seesaw 

In the supersymmetric inverse seesaw (ISS) [148-150], the MSSM particle content is ex¬ 
tended with 3 generations of right-handed neutrino superhelds and 3 generations of 
singlet superhelds X The superpotential takes the form 

W = Wmssm + , (16) 

where we have omitted family indices. W and are general 3x3 complex mass matrices 
and /Xi/ is a complex symmetric 3x3 matrix. While M/j generates a lepton number 
conserving Dirac mass term for the fermion singlets, violates lepton number by two 
units. This Majorana mass term also leads to a small mass splitting in the heavy neutrino 
sector, which is then composed by three quasi-Dirac neutrinos. The corresponding soft 
SUSY breaking Lagrangian is given by 

+ {T^N^LH^ + Bm^N-X + ^B^^XX + X^m^-N^ + h.c.), (17) 

where Bmj^ and B^^ are the new parameters involving the scalar superpartners of the 
singlet neutrino states. Notice that while the former conserves lepton number, the latter 
violates lepton number by two units. Finally, contains the soft SUSY breaking 

terms of the MSSM. 

The scalar potential of the model is such that the neutral components of the Higgs 
superhelds get non-zero VEVs, 

(//s> = Y . (■»"> = Y. ( 18 ) 


®More minimal realizations of the ISS are possible [151-155]. However, for simplicity, we will stick to 
the most common version with 3-1-3 singlet superhelds. 
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triggering electroweak symmetry breaking (EWSB). This induces mixings in the neutrino 
sector. In the basis v = , X), the 9x9 neutrino mass matrix is given by 

/ 0 0 \ 

-Adiss = j 0 Mji j . (19) 

\ 0 j 

where mD = -^YyVu- Assuming the hierarchy /i^ mo Mr, the mass matrix Adiss 
can be approximately block-diagonalized to give the effective mass matrix for the light 
neutrinos [156] 

’T^iight - rn^ /ij, mn • (20) 

On the other hand, the other neutrino states form three heavy quasi-Dirac pairs, with 
masses corresponding approximately to the entries of Mr. 

Eq. (20) has the same form as Eq. (15), with Mss ~ Mr. Therefore, by taking a small 
parameter, the model allows for small neutrino masses, sizable Y^, Yukawa couplings 
and singlet neutrinos at the TeV scale (or below). 

4.2 LFV in the supersymmetric inverse seesaw 

The presence of light singlet neutrinos induces all sorts of effects. Here we will concentrate 
on their contributions to LFV processes. For some recent works on phenomenological 
aspects of light singlet neutrinos see [120,123,126,130,131,157-165]. 

There is a vast literature on LFV in models with light singlet neutrinos. Potentially 
large enhancements, with respect to the usual high-scale models, were already pointed out 
in early studies [55,69,75,150]. More recently, several works have explored in detail the 
LFV anatomy of these models, highlighting the relevance of (non-SUSY) box diagrams 
induced by singlet neutrinos [145,166-168], computing Higgs penguin contributions [169] 
and showing enhancements in the usual photon penguin contributions [159]. Regarding 
purely supersymmetric contributions, the relevance of Z penguins with right sneutrinos was 
recently readdressed in [170], solving an inconsistency in the analytical results of [57] and 
[171]. Finally, [172] constitutes the hrst complete analysis of LFV in the supersymmetric 
inverse seesaw, taking into account all possible contributions, supersymmetric as well as 
non-supersymmetric. 

We will now present the main results in [172]. These were obtained using FlavorKit 
[173], a tool that combines the analytical power of SARAH [174-178] with the numerical 
routines of SPheno [179,180] to obtain predictions in a wide range of models, based on 
the automatic computation of the lepton flavor violating observables. See [181] for a 
comprehensive and pedagogical review of this set of tools. 

In the following, we will discuss numerical results obtained using universal boundary 
conditions at the gauge coupling unihcation scale, mcuT — 2 ■ 10^® GeV, setting Mrusy = 
mo = Mi /2 = — Aq. In addition, we hxed = 100 /ij,, Bm^ = ^00 Mr, tan/3 = 10, /i > 0 
and considered a degenerate singlet spectrum (M]j = Mr with i = 1,2,3). Furthermore, 
we hxed the V, Yukawa couplings using a modihed Casas-Ibarra parameterization [86], 
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Figure 4: BR(/i —)■ ej) as a function of Msusy and Mr. On the left-hand side Msusy = 1 
TeV is hxed, whereas on the right-and side we set M^i = 2 TeV. The other parameters are 
given in the text. The gray area roughly corresponds to the parameter space excluded by 
the LHC SUSY searches. Figure taken from [172]. 



adapted for the inverse seesaw [157,182], in order to reproduce the neutrino squared mass 
differences and mixing angles observed in oscillation experiments [183] (see also [7] for an 
update). 

A general conclusion one can draw from [172] is that the LFV phenomenology strongly 
depends on Mr and Msusy- The hrst scale determines the mass of the singlet neutrinos, 
whereas the second one sets the superparticle masses and their relative size determines the 
phenomenology. This can be seen in Fig. 4, where BR(/i —)■ ey) is shown as a function 
of Msusy and Mji. The results are displayed in three curves: the full observable, the 
SUSY contributions and the non-SUSY ones. The latter consist of contributions from u- 
and u-H^ loop diagrams, thus involving the singlet neutrinos in combination with the 
W boson or a charged Higgs. One hnds that the relative weight of SUSY and non-SUSY 
contributions is given by the hierarchy between these two mass scales. For Msusy S> M^, 
non-SUSY contributions induced by the singlet neutrinos dominate the /i —)■ ey amplitude, 
whereas for Msusy ^ Mr, the usual MSSM contributions generated by chargino/sneutrino 
and neutralino/slepton loops turn out to be dominant. Moreover, we hnd that non-SUSY 
contributions can have strong cancellations. 

A similar behavior is found in case of the 3-body decays ii —)■ 3ij. In Fig. 5 we display 
numerical results for BR(/i —>■ 3e) as well as for various contributions to this observable. 
The anatomy of this decay is more involved, since more types of Feynman diagrams con¬ 
tribute to the amplitude: SUSY and non-SUSY photon, Z and Higgs penguins, as well 
as box diagrams. As for U —)■ ij'j, we observe that non-SUSY contributions dominate for 
low Mr (< Msusy)- In particular, we see on the left-hand side of Fig. 5 that non-SUSY 
boxes become completely dominant as soon as one goes to Mr values below ~ 2 TeV. 
This generic feature was already noted in [145,166-168]. For higher values of Mr SUSY 
contributions, and in particular the standard photon dipole penguin, dominate. On the 
right-hand side we hnd complementary information, with the different contributions as a 
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Figure 5: BR(/i —3 e) as a function of Msusy and M/j. On the left-hand side Msusy = 1 
TeV is fixed, whereas on the right-and side we set M^i = 2 TeV. The other parameters are 
given in the text. The gray area roughly corresponds to the parameter space excluded by 
the LHC SUSY searches. Figure taken from [172]. 


function of Msusy for a fixed M/j. It is worth noticing that supersymmetric Z penguins 
never dominate. 

Finally, analogous results are obtained for the /i — e conversion in nuclei rates. Inter¬ 
estingly, the large non-SUSY boxes found at low Mr break the dipole dominance, leading 
to a clear departure from Eq. (3). This is illustrated in Fig. 6, where BR(/r —)■ ey), 
BR(/i —)■ 3 e) and the /i — e conversion rates in Ti and A1 are shown as a function of Mj^. 
Indeed, for Mr < 500 GeV, the rates for all LFV processes have similar sizes. In this 
scenario, experiments looking for /i —)■ 3 e and /i — e conversion in nuclei will soon provide 
the most stringent constraints in this model. 


5 R-parity violating models 


The particle content and symmetries of the MSSM allow for the following superpotential 
terms 




2^^ijkLiLjCj^ -|- X—j^L{Qjdf^ -|- SiLiHy^ 


( 21 ) 


where we have explicitly introduced family indices. Here and are the right-handed up- 
quark and charged lepton superfields, respectively, and the rest of the superfields have been 
already dehned. The hrst three terms in break lepton number (L) whereas the last 
one breaks baryon number (B). In principle, these couplings are not welcome, since they 
give rise to many lepton and baryon number violating processes, never observed in nature. 
For example, the simultaneous presence of the X' and X" couplings would lead to proton 
decay [184,185]. This phenomenological problem is solved in the MSSM by introducing by 
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Figure 6 : BR(/i —)■ 67 ), BR(/i —)■ 3 e) and /i — e conversion rates in Ti and A1 as a function 
of Mr. Msusy is fixed to 1 TeV. The gray area roughly corresponds to the parameter 
space excluded by the LHC SUSY searches. Figure taken from [172], 


hand a new discrete symmetry that forbids all terms in hF^p . This symmetry is known as 
R-parity [3,4] and is dehned as 

Rp = ('_1)3(S-L)+2 s _ ^ 22 ) 

Here s is the spin of the particle. It is straightforward to verify that all terms in (21) 
break R-parity and thus they are forbidden once this symmetry is imposed. The MSSM is 
dehned as R-parity conserving. 

However, several arguments can be raised against R-parity: 

• R-parity is imposed by hand. Unlike the SM, where L and B conservation is 
automatic, in the MSSM this has to be forced by introducing a new symmetry, not 
derived from hrst principles. This is clearly a step back from the SM. 

• R-parity does not solve fast proton decay. It is well known that R-parity does 
not forbid some dangerous dimension-5 operators that lead to proton decay [186-188]. 
For example, the operator O 5 = -^QQQL has Rp{0^) = -|-1 and thus conserves 
R-parity. The bounds obtained from the non-observation of proton decay imply 
that, even for M = Mpianck, / must be smaller than 10“^ [189]. In order to forbid 

and other similar dimension-5 operators, one may resort to additional havor 
symmetries [190]. 

• There is no reason to forbid all the L and B violating operators. Proton 
decay requires the simultaneous pesence of L and B violating couplings. Therefore, 
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it is sufficient to impose the conservation of just one of these two symmetries in 
order to forbid proton decay. This has led to the consideration of alternative discrete 
symmetries which allow for either L or B violation while protecting the proton. An 
example of such symmetries is baryon triality {Z^) [186,191]. 

Furthermore, there are several good motivations to consider R-parity violating {fip) 
scenarios. The violation of lepton number by any of the hrst three couplings in Eq. (21) 
automatically leads to non-zero neutrino masses [192-194]. Moreover, the presence of 
couplings leads to a rich collider phenomenology due to the decay of the LSP. This can be 
translated into longer decay chains, changing the expected signatures at the LHC [195,196]. 
In fact, fi,p has also been considered as a way to relax the stringent bounds on the squark 
and gluino masses, see for example [5,6,197,198]. 

Finally, in l^pthe standard neutralino LSP is lost as a dark matter candidate. Therefore, 
alternative candidates must be considered. Examples in the literature include (i) gravitinos 
[199-201], (ii) the axion [202,203] or (hi) its superpartner, the axino [204,205]. For general 
reviews on R-parity violation and collections on bounds on the couplings see [206-209]. 

We will now discuss separately the LFV phenomenology of two very different supersym¬ 
metric scenarios with R-parity violation: explicit R-parity violation {e-fi,p) and spontaneous 
R-parity violation (s-fip). 

5.1 Explicit R-parity violation 

The most characteristic signatures of R-parity violating models are, of course, processes 
with L or B violation. Nevertheless, processes that violate lepton flavor can provide inter¬ 
esting signatures as well and, in fact, they can be more attractive due to the large number 
of upcoming LFV experiments. 

Higgs LFV decays 

After the historical discovery of the Higgs boson [1,2], a lot of effort has been put into 
the determination of its properties. In particular, the Higgs boson decays may contain a 
lot of valuable information, with potential indications of new physics. Recently, the CMS 
collaboration reported on an intriguing 2.4cr excess in the h —)■ r/x channel [32]. This hint, 
which translates into BR(h —)■ r/x) = (0.84 ^q 37), is based on the analysis of the 2012 
dataset, taken at y/s = 8 TeV and an integrated luminosity of 19.7 fb“^. This large Higgs 
LFV ® branching ratio is quite challenging and most NP models cannot accommodate 
it [212]. In fact, the flavor conserving Higgs decay h —)■ rr has a branching ratio of 
only ~ 6%, not much higher than the LFV one found by CMS. Although independent 
con&mation by ATLAS, as well as additional statistics in CMS, would be required in 
order to promote this hint to the category of evidence of new physics, it is interesting to 
explore different models in order to determine what type of frameworks can accommodate 
this signal. 

®For pioneer works on Higgs LFV decays see [210,211]. 
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Figure 7: Tree-level contributions to —)■ r/r. On the left, a. B X contribution. On the 
right, a e] (denoted as fij in [213]) contribution. Figure borrowed from [213]. 
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Figure 8 : BR(h —)■ fir) contours on the B 2 — A 232 plane. The continuous horizontal and 
vertical lines show approximate limits due to neutrino masses (in case of B 2 ) and charged 
current experiments (in case of A 232 )- Figure borrowed from [213]. 
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Regarding supersymmetric models, several scenarios have been recently explored, some 
of them even before the CMS hint was announced. In particular, the authors of [213,214] 
considered an extension of the MSSM including all L violating couplings in Eq. (21). The 
particles-sparticles mixing due to the couplings induce Higgs LEV decays at tree-level, 
thus potentially being able to reach branching ratios as high as the one found by the 
CMS collaboration. Two specihc examples are shown in Fig. 7. However, the existing 
experimental bounds on the relevant combinations of couplings contributing to h —)■ 
rp forbid such large LEV branching ratios. For example, in the B\ contribution both 

parameters are strongly constrained. In case of Bi, the ftp mixings between the Higgs 
boson and the sneutrinos {C D BiLiHy), they have strong bounds since they induce non¬ 
zero neutrino masses [215-217]. Moreover, the A couplings are constrained by charged 
current experiments [206]. Once these constraints are taken into account, the maximum 
BR(h —)■ Tfi) one can get is not very impressive. This is illustrated in Fig. 8, where 
BR(h — )■ r/i) contours are drawn on the B2 — X232 plane. From this hgure one concludes 
that BR(h —)• rp) can reach, at most, a few x 10“®, clearly below the CMS hint. Similar 
conclusions are obtained when other combinations of ftp parameters are considered. 

The supersymmetric inverse seesaw has also been considered as a possible setup to 
reproduce a Higgs LFV branching ratio into r/i at the 1% level [218]. In this case, h —)■ r/i 
takes place at 1-loop, naturally suppressing the branching ratio. As a result of this, as well 
as due to the constraints from other LFV processes such as it —>■ ij'y, one finds that the 
maximum allowed BR(h. —)■ Tfi) is about ~ 10“®. Therefore, this model cannot account 
for a branching ratio as obtained by CMS either. In order to conclude this discussion 
on h —)■ r/i with a positive note, let us mention that known models that can account for 
BR(h —>■ r/i) ~ 1% exist in the literature. They all involve extended Higgs sectors. In 
particular, it has been shown that Two-Higgs-Doublet models of type-HI, in which both 
Higgs doublets can couple to up- and down-type fermions, can easily accommodate the 
CMS signal [212,219-223]. The MSSM, being a Two-Higgs-Doublet models of type-H, 
cannot 

Trilinear R-parity violation and LFV 

In principle, the usual LFV processes studied in R-parity conserving models can be studied 
in the R-parity violating ones and, in some cases, they get additional contributions. This 
is the case of £* — )■ 3 ij, M —)■ ii£j and r — )■ M£i, where M is a neutral meson, which, in the 
presence of trilinear ftp couplings, can be induced at tree-level [227]. This is represented 
in Fig. 9, where two examples are shown: £i —)■ 3£j induced by A couplings and sneutrino 
exchange and r —)■ M/i induced by A' couplings and squark exchange. This allows one to 
derive a large collection of bounds on the size of the trilinear couplings and the masses of 
the superparticles mediating the LFV decays [228,229]. 

Let us consider the right-hand side of Fig. 9. After dressing the quarks in the hnal 

fact, several studies have shown that one cannot accommodate BR(/i —)■ r/i) ^ 1% when the 
low-energy theory is the MSSM [224,225]. The same conclusion applies to heavy Higgs LFV decays [226]. 
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Figure 9: To the left, li —)• Sf'j induced by trilinear A couplings and sneutrino exchange. 
To the right, t ^ Mfi induced by A' couplings and squark exchange. 

state, this Feynman diagram induces r —)■ Mfi at tree-level. Since the exchanged particle, 
a down-type squark in this case, is much heavier than the rest of particles, this process can 
be well described by the 4-fermion effective Lagrangian 



(23) 


obtained after integrating out the heavy squark. One can now use this Lagrangian and, 
together with the relevant hadronic form factors, compute rates for processes such as 
T —)■ /i7r+7r“. The authors of [229] followed this method and used Belle results on searches 
for r LFV decays [230-232] to obtain the limit 



(24) 


One can exploit this idea using other LFV observables involving mesons. We refer 
to [228,229] for a more complete list of constraints. 

Similarly, trilinear couplings can also trigger /i — e conversion in nuclei, induced by 
diagrams very similar to the one on the right-hand side of Fig. 9. Interestingly, in this case 
fi — e conversion in nuclei would take place at tree-level, while the more popular /i —)■ ey 
would take place at 1-loop. This has been recently pointed out by the authors of [233], 
who argue that experiments looking for /i — e conversion in nuclei might be the hrst (and 
perhaps the only ones) to observe a non-zero signal in the next round of experiments. 

Other results on LFV in fip scenarios 

Before concluding, let us briefly comment on other aspects of LFV in ^p models. An 
interesting feature of ^p models is that some lepton number violating processes at colliders 
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Figure 10: Sneutrino decay to ikij in bilinear R-parity violation. The open circle with 
a cross indicates the fip induced mixing between charginos and charged leptons. Figure 
taken from [234]. 


might look like lepton flavor violating ones. This is for example the case of sneutrino decay 
in bilinear R-parity violation [234], as shown in Fig. 10. This process is possible thanks 
to the mixing between the MSSM charginos and the standard charged leptons. At the 
LHC, if the sneutrinos are directly produced, the absence of missing energy would make 
this process look like a LFV one. 

Finally, for other recent works on LFV in fi,p models see [235-237]. 

5.2 Spontaneous R-parity violation 

A very attractive scenario for LFV is that of spontaneous R-parity violation. When a scalar 
held oddly charged under R-parity gets a VEV in a theory with a R-parity conserving 
Lagrangian, R-parity gets spontaneous broken. Here we will concentrate on spontaneous 
L and Rp violation. Although R-parity is a discrete symmetry, its breaking comes along 
with the breaking of the continuous global symmetry U{1)l- This implies the existence of 
a massless Goldstone boson, usually called the majoron (J) [238,239]. 

The nature of the majoron is crucial for the phenomenological success of the model. 
In fact, in the hrst model with s-j7p [240], the breaking of R-parity was triggered by 
the VEV of a left-handed sneutrino. This simple setup was eventually excluded since the 
doublet nature of the majoron leads to conhict with LEP bounds on the Z boson invisible 
decay width and astrophysical data [203,241]. However, more rehned models where the 
violation of lepton number is induced by a gauge singlet are perfectly valid possibilities. 
As a benchmark example of this family we will consider here the model introduced in [242]. 
For alternative models with gauged lepton number see for example [243-245]. 

In the model of Ref. [242] , the particle content is extended with three additional singlet 
superhelds, namely, S and <h, with lepton number assignments of L = —1,1,0 respec¬ 
tively. By assumption, the Lagrangian of the theory conserves lepton number. Therefore, 
the superpotential can be written as 

WsRPv = VFmssm + - hoH^Hj + + ^5^ . (25) 
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For simplicity, one can simply consider one generation of and S superfields. Several 
scalar fields acquire VEVs after electroweak symmetry breaking. In addition to the usual 
MSSM Higgs boson VEVs, Vd and these are {$) = Vis,j\/2, {N^) = (S) = vsl\/2 

and = Vi/\/2. This vacuum conhguration breaks lepton number and R-parity. In fact, 
we notice that vr ^ 0 generates the effective bilinear j^^terms e* = Y/vr/ \/2. Furthermore, 
neglecting Vi -C Vr, Vr, one hnds the resulting majoron prohle 

, ( 26 ) 


where V = + n|. Eq. (26) shows that the majoron inherits the singlet nature of the 

scalar helds that break lepton number with their VEVs, thus suppressing the couplings to 
the Z boson and evading the stringent LEP bound. 

Here we are interested in novel LFV features due to the presence of the majoron 
This new massless state dramatically changes the phenomenology both at collider and low- 
energy experiments [248,249]. In particular, it leads to new LFV processes, such as p —)■ eJ 
or p —)■ eJy. The exotic muon decay /i —>■ eJ was hrst studied in [250] and later revisited 
in [249], where the decay with an additional photon was also considered. Furthermore, the 
impact of the majoron on /i — e conversion in nuclei was discussed in [251] 

The rate of the /i —)■ e J decay is determined by the e — p — J coupling, Oe^j, which, in 
the model under consideration, is of the form Oe^ij ^ ^ ^ RPV parameters. This makes 
us conclude that, in general, one expects large partial muon decay widths to majorons for 
low Vr. However, currently there are no experiments looking for /i —)■ eJ and the current 
best limit on the branching ratio, BR(p —>■ eJ) < 10“®, dates back to 1986 [253]. Regarding 
the decay including a photon, /x —>■ eJy, one can proht from the MEG experiment and its 
search for the more popular channel /x —>■ ey. 

The two branching ratios are related by 


cx 

BR(^ t eTy) — ymin) BR(p, y eT). 

zn 

Here I{xmin,ymin) is a 3-body phase space integral dehned as 


(27) 


Yi^XYnim ymin) I dxdy 


{x - 1)(2 -xy-y) 
x/2(l -x-y) 


the dimensionless parameters x, y are dehned as 


2E, 


X = 


m„ 


y = 


2R 


m, 


(28) 


(29) 


and Xmin and ymin are the minimal electron and photon energies that a given experiment 
can measure. Indeed, the integral in Eq. (28), which would contain infrarred and collinear 
divergences, is regularized by the specihc choices made by an experiment. 

^^Another interesting signature present in majoron models is the invisible decay of the Higgs boson, 
JJ [246,247]. 

^^See also [252] for similar LFV processes in the context of invisible axions. 
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Figure 11: The phase space integral for the decay /r —>■ eJy as a function of Xmin for three 
different values of ymin = 0.95, 0.99, 0.995 from top to bottom and for two different values 
of cos^e..;.. To the left cos9e-y = —0.99, whereas to the right cos 6^67 = —0.99997. Figure 
taken from [249]. 


As explained above, the main advantage of p —)■ eJy is the existence of the MEG 
experiment. However, the question is whether it is sensitive to this exotic LFV process or 
not. Fig. 11 shows the value of the phase space integral X{xmin, ymin) as a function of Xmin 
for three different values of ymin and for two choices of cos 6*67 (the relative angle between 
the electron and photon directions). Unfortunately, the MEG experiment is specihcally 
designed for a single search. In fact, the cuts used in the search for /i —>■ ey are very 
restrictive: Xmin > 0.995, ymin > 0.99 and [vr — 9e^\ < 8.4 mrad. For these exact values 
one hnds a tiny phase space integral, X ~ 6 ■ 10 “^°. As a consequence of this, a limit 
for BR(/i —)■ ey) of the order of < 10“^^ would translate into the useless limit BR(p —)■ 
eJ) < 0.14. To improve upon this bound, it is necessary to relax the cuts. For example, by 
relaxing the cut on the opening angle to cos 6^67 = —0.99. However, this is prone to induce 
additional unwanted background events. In particular, accidental background from muon 
annihilation in flight. Therefore, although one could in principle increase the value of the 
phase space integral X{xmin,ymin), the background in that case would make the search for 
a positive signal impossible. This discussion suggests that a better timing resolution of the 
experiment would be welcome in order to reduce the background and be sensitive to hnal 
states including majorons. 

6 Summary and conclusions 

In summary, we have reviewed the lepton flavor violating phenomenology of several non- 
minimal supersymmetric models: high-scale and low-scale seesaw models as well as models 
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with explicit or spontaneous R-parity violation. The main conclusion from this overview 
is that the lepton flavor violating signatures can be very different from those found in the 
MSSM. This translates into two important messages: 

• For the theorists: Lepton flavor violation might be much more intricate than what 
minimal models predict. Therefore, we should be careful when extrapolating our 
expectations (derived from the MSSM) to extended frameworks. 

• For the experimentalists: Although minimal models are of course well motivated, 
lepton flavor violation might show up in non-standard channels. We must be ready 
to avoid missing a relevant signal. 

Properly identifying the underlying physics will be crucial in case a positive observation 
in one or several LFV experiments is made. This problem might be soon have to be 
addressed, given the exciting projects that are currently going on or soon starting their 
search for LFV. Hopefully, this review, as well as the many phenomenological studies in 
the bibliography, will help shedding some light on this matter. 
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